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Abstract 

A  study  has  been  made  of  the  feasibility  of  accurate  numerical 
determinations  of  the  transmission  of  gamma  rays  through  large  thicknesses 
of  materials.  The  first  procedure  investigated  consists  in  regarding  the 
total  probability  of  photon  transmission,  N^,  as  the  sum  of  the  probabili¬ 
ties  Nn,  where  Nn  is  the  probability  of  photon  transmission  with  exactly  n 
scatterings.  The  total  expected  transnitted  energy,  E^,  is  similarly 

considered  to  be  given  byZjE  .  A  numerical  calculation  of  N  and  E„  has 

n  n  n 

been  made  for  n  «  0,  1,  2,  3  for  a  slab  of  uranium  20cm  thick  upon  which 

2 

photons  are  incident  normally  with  energy  o  ■  10  me  •  The  maximum  value 

of  n/N  occurs  at  n  «  2  and  of  E  /E^  at  n  -  1.  These  calculations  are 
rr  0  n  o 

also  adapted  to  a  slab  of  lead  35cm  thick.  Consideration  has  been  given 

to  the  behavior  of  N  and  E  for  large  n,  and  estimates  are  thereby  made 

n  n 

for  N,  and  E.  •  The  second  procedure  consists  in  deriving  the  transmission 

V  V 

through  a  thick  slab  from  a  succession  of  transmissions  through  thin  slabs. 
The  transformation  of  an  incident  photon  distribution  into  the  distribution 
transmitted  through  a  thin  slab  is  conveniently  expressed  as  a  matrix,  and 
the  total  transmission  is  then  given  by  the  iteration  of  the  matrix  on  the 

successive  transmitted  distributions.  Numerical  results  obtained  by  this 
procedure  for  particular  incident  photon  distributions  are  presented. 


I.  INTRODUCTION 
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The  purpose  of  the  present  investigation  has  been  a  study  of  the 
feasibility  of  straightforward  numerical  determinations  of  the  trans¬ 
mission  of  gamma  rays  through  large  thicknesses  of  materials  containing 
heavy  elements.  The  limitation  of  the  discussion  to  heavy  elements  makes 
possible  sane  simplifications  essentially  because  of  the  large  probability 
of  absorption  by  the  photoelectric  effect  of  photons  which  are  degraded 
to  low  energy.  The  present  approach  to  the  gamma  ray  transmission  problem 
is  based  on  the  notion  that  one  need  consider  only  those  transmitted  gamma 
rays  which  have  suffered  relatively  few  scatterings  even  fbr  thicknesses 
of  materials  of  approximately  20  mean  free  paths.  The  validity  of  this  view 
is  demonstrated  by  the  results  which  are  presented  in  the  following  sections. 
Specific  calculations  have  been  made  of  the  attenuation  of  ganma  rays  with 
incident  energy  10  me  through  thicknesses  of  uranium  up  to  20cm.  These 
calculations  are  adjusted  to  give  the  attenuation  through  thicknesses  of 
lead  up  to  35  cm. 

The  elementary  processes  of  gamma  ray  interaction  with  matter  that 
are  taken  into  account  are  the  photoelectric  effect.  Compton  scattering, 
and  pair-production.  Sinoe  the  numerical  results  for  gamma  ray  trans¬ 
mission  will  depend  upon  the  values  taken  for  the  absorption  coefficients 
for  these  processes,  a  partial  table  of  the  values  for  uranium  and  lead 
used  in  the  present  calculations  is  given  in  Table  I.  If  a  minor  alteration 
is  made  in  these  values,  the  effect  on  the  transmission  values  given 
here  could  be  determined  without  great  difficulty.  The  probability  of 
Compton  scattering  is  assumed  to  be  given  by  the  Klein-Nishina  formula 
in  which  the  effects  of  polarization  of  the  radiation  have  been  averaged 
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out.  Thus,  the  partial  polarisation  which  arises  upon  scattering  of 
originally  unpolarised  radiation,  and  the  alteration  in  the  probability 
of  further  Compton  scattering,  has  been  disregarded.  The  gamma  ray  energies 
of  present  concern  are  sufficiently  high  so  that  these  polarisation  effects 
are  believed  to  be  unimportant.  The  possible  contribution  to  the  trans¬ 
mitted  radiation  from  nBremssthralungn  produced  by  Compton  recoil  electrons 
is  not  included  in  the  present  calculations]  the  range  of  gamma  ray  energies 
considered  here  is  known  to  be  too  low  for  this  contribution  to  be  signifi¬ 
cant. 

Two  methods  for  the  determination  of  the  probability  of  transmission 
of  a  photon  through  a  slab  of  material  of  thickness  a  will  be  considered  here. 
In  the  first  method,  the  total  probability  of  transmission  is  taken  as  the 
sum  of  the  probability  of  transmission  with  no  scattering,  plus  the  proba¬ 
bility  of  transmission  with  one  scattering,  plus  the  probability  of  trans¬ 
mission  with  two  scatterings,  etc.  In  the  second  method  the  slab  of  thickness 
Jl  is  divided  into  a  series  of  thin  slabs  and  the  transmission  through  the 
total  thickness  is  determined  from  the  transmission  through  the  series  of  thin 
slabs  in  succession. 

II  •  Successive  Scattered  Contributions  to  the  Transmitted  Radiation 

Consider  a  homogeneous  slab  of  material  which  has  infinite  extent  in  the 
y  -  and  z  —  directions  and  which  has  thickness  a  in  the  x  —  direction, 

0  <  x  <  a.  One  may  readily  give  an  integral  expression  for  the 
probability  that  a  gamma  ray  incident  on  the  face  x  -  0  is  not 


-  3  - 
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absorbed  in  the  slab  and  emerges  after  exactly  n  collisions  from  the  face 

x  ■  a.  Let  the  gamma  ray  enter  with  the  energy  ao  in  a  direction  making 

an  angle  «J»Q  with  the  normal  to  the  face  of  the  slab,  travel  a  distance  sQ 

within  the  slab  to  the  first  collision,  scatter  (if  it  is  not  absorbed) 

into  a  direction  making  an  angle  ©1  with  the  original  direction  and  an 

ar*fil®  with  the  normal  to  the  slab  face,  travel  a  distance  s^,  to  a 

second  collision,  scatter  into  a  direction  specified  by  ©2  and  <^,  and  so 

on  (cf.  Fig.  1).  After  n  collisions,  the  path  angles  are  0  and  <1.  ,  and 

'  n  n' 

■n  is  the  distance  traveled  from  the  (n-l)th  collision  to  the  exit  face  of 
the  slab  at  x  ■  a,  The  path  line  after  the  k'th  collision  is  characterised 
by  an  azimuthal  angle  0k  in  addition  to  the  angles  ©k  and  these  angles 
are  connected  by  the  familiar  relation 

COS  *k+l  “  CO8^kcos0k+l  *  8in  +k  8in\+l  C08^k+lJ  k*0*  1»  z»  •  •  • » ( n-1 ) . 
The  energy  of  the  photon  between  the  k'th  and  the  (k+l)th  collision  is  ak 
in  units  of  me  ,  and  the  value  of  the  total  absorption  coefficient  for 
this  energy  ak  will  be  denoted  by  The  successive  values  of  the  energy 
ak  are  related  by  the  Compton  formula 1 

„  gk-l  . 

k  1  -  a^C  1-cos  «k) 

The  probability,  Nn,  that  the  photon  will  be  transmitted  after  exactly 
n  collisions  by  any  possible  path  is 
o  n 

Vr2  ' 

n  1  2 


N 


-li  s  n-1 
ti  n 


•  •  • 

c 


k-0 


~*1k8k  ,  w  ,Q 

e  Vi  d3k  dsk*i, 


0) 


where  Vis  the  number  of  electrons  per  enr  in  the  slab  material,  r  is  the  electron 


radius,  e  /me  ,  and  qk+^  is  obtained  from  the  Klein-Nishina  differential 


cross-section  as: 


-  4  - 
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q, 


lk+l 


sin  ©ktl 


i-  =os2  e.  .  .  “k  U-co»  e^)2 
r*\Ti^o.  s^j.) 


l*ak(l-co8  ©k+1); 

The  3n-dimensional  apace  S  over  which  the  integration  in  Eq.  (1) 
must  be  performed  needs  detailed  description.  A  point  (so,©1,  ^1 
“n-l1  ®n*  roPre8ents  a  path  followed  by  the  photon  in  reaching  the  n'th 
collision  and  is  a  possible  path  or  a  point  in  S  if  all  these  collisions  are 
within,  or  at  the  boundary  of,  the  slab.  Any  value  of  0^  and  ©k  is  possible 
as  long  as 


0  <K  <  2*  » 


0  <  ek  <  n  . 


k  ■  1,  2,  •  .  «,n« 


The  first  collision,  however,  will  not  be  in  the  slab  unless 

0  <  s  cos  <i  <  a. 

—  o  o  — 

Further,  the  value  of  the  angle  determines  whether  the  first 
. collision  has  given  forward  or  backward  scattering.  If  the  scattering  is 
forward, 


0  <  ^  <  */2  , 

then  the  second  collision  will  be  possible  only  if 

0  <  s^  cos  <|^  <  a-sQ  cos  . 

If  the  scattering  is  backward, 

n/2  <  4^  <  Tt  , 

then  the  inequality 

0  >  s.  cos  <U  >  -  s  cos  4 
—  1  1  —  o  To 

must  hold.  Continuing  in  this  way,  one  sees  that  the  space  S  is  defined 
by  the  inequalities j 

0  <0k  <  2n  ;  °  <  ©k  <  ti  j  (2a) 
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k-2 

0  -  *k-l  008  \-l  -  a  -  EjJ  cos  4^  , 

if  °s  Vis! > 

(2b) 

k-2 

°2Vi  c"  <W  2  "  *i  c°*  «i  * 

if  2  <  V-l  -  "  • 

(2c) 

For  each  Inequality  of  (2),  k  ■  1,  2,  ...,  nj  it  is  to  be  noted  that  for  k  »  1, 
k-2 

one  replaces  s^  cos  4^  by  zero,  the  last  inequality  (2c)  which  gives  the 
range  of  s^^  cos  4’j^  is  trivial,  and  the  inequality  defining  the  range  of 
4^2  loses  significance. 

The  space  3  divides  naturally  into  2n”^  subspaces,  where  each  subspace 
represents  a  particular  sequence  of  forward  and  backward  scatterings.  The 
scatter  forward  or  backward  on  the  last,  or  n'th,  collision  does  not  figure 
in  these  inequalities,  but  does  enter  by  way  of  the  integrand.  For  transmission, 
sn  is  defined  by  the  relation 


8n  co»  «  a  -  ZI  Si  cos  ^  ,  if  0  <  *n  <  |  , 


and  for  reflection  by 

n-1 

«n  008  *n  "  ”  ZZ  s±  cos  ^  ,  if  |  <  ^  <  7t  .  (3b) 

i»0 

Hence,  the  number  of  subspaces  is  increased  to  211,  and  one  may  expect  that  the 
calculation  of  Nn  as  a  practical  matter  will  require  separate  treatment  for 
each  sub space. 


The  integral  formula  for  Nn  does  not  appear  tractable  to  analytic  treatment 

of  the  integration  over  the  angle  variables,  but  it  will  be  noted  that  the 

integration  with  respect  to  the  s^'s  may  be  readily  performed.  In  order  not  to 

single  out  a  particular  subspace,  one  may  combine  Eqs.  (3a)  and  (3b)  into  the 

one  relation  , 

n-1 

Bn  cos  “  an  “  g,  si  cos  * 

■where  aQ  is  assigned  the  value  a  or  0  according  as  (3a)  or  (3b)  applies. 
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Similarly,  the  integration  limits  on  the  variable  s,  as  given  by  (2b)  or 

k-2  K 

(2c)  will  be  taken  to  be  0  and  0o®  <k^,  where  a^  has  the  value 

a  or  0  according  as  the  subspace  defined  by  (2b)  or  (2c)  applies.  Then 
Eq.  (1)  requires  the  integration  of  the  function 


“P  -  £  ’1  -  -  **[-  £  H  ‘ ij 

n-2 

first  with  respect  to  8n-1  over  the  interval  (0,  a^  -  g  8l  cos  4^),  then 

rW? 

with  respect  to  a  2  over  the  interval  (0,  a  ,  -  >  8<  008  <ki )»  etc.  This 

i»0  1  1 

series  of  operations  has  a  recursive  nature  which  can  be  made  evident  by  the 
following  device.  Let 

f  T  k  7  r  k  1 

exp  -  k.(a.  -  2_J  •  i  cos  O/  •  exp  -  Y2  p s,  /  for  k  >  0; 
LJJi«01  XJ  /  ioQ  1  1  I3 


BprVj  h 


for  k  *■  -1 . 


Then,  if  represents  the  operation  of  integrating  with  respect  to  s.  between 

k-1  k 

the  limits  0  and  a.  -  >  s.  cos  4.  ,  one  has 

i-0  1  1 


[F  -e  ;  V  ' 

\  J’Z-f  /s  ,  \ 

(  /y  'Ay  y  Ctfii 


Further  iterations  proceed  similarly.  It  is  evident  that  the  s-integrations 
in  Eq.  (1)  are  the  end  of  a  sequence  of  these  operations  of  which 


-  7  - 
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(4a)  and  (4b)  are  the  first  two;  the  sequence  in  is  obtained  from  (4) 

by  these  substitutions:  k  =  n-1,  1  =  n,  X  =  u  /cos  4  .  For  n  =  1,  one  has 

j  >  o  j  n  nr  Tn  * 


-X,  a.,  -X,  (a, -a)  -X  a 

11  11  o 

e  -e  e 


f  t - 7 -  q,  dO,  d0,  , 

X  -X, )  cos  (I  1  1  HL  ’ 

0  1'  To 


and  for  n  =  2 


N2  ^al,a2^  “I'T' 


rr  e^VX^’-a)=-X< 


,Xr4a)  cos  4c 


te-Ma;-aj)  i 


-X,af  -X , ( a,  -a )  -X„a 


I X 0  —X 1 )  cos 


^1  “  C0S  ^1 


1lq2d^ldeid^2d°2  * 


Thus,  the  probability  of  a  photon  passing  through  the  slab  with  exactly  one 
collision  is  given  by  Eq.  (5)  as  N^(a),  and  the  probability  of  being  reflected 
out  through  the  incident  face  of  the  slab  with  one  collision  is  N^(0).  The 
probability  of  transmission  with  exactly  two  collisions  is  given  by  Eq.  (6) 
as  the  sum  of  the  two  probabilities  ^(a,  a)  and  ^(0,  a),  and  the  proba¬ 
bility  of  reflection  out  through  the  incident  face  with  two  collisions  is 

the  sum  of  the  two  probabilities  N^a,  0)  and  ^(0,  0).  If  one  defines 
Vi  =  *ia>  bben  these  probabilities  may  be  written  as  follows: 


i  2  -v 
\  v  r  a  o 
Ni(a)  =  — —  e 


N,(0)  = 


ff  "(vi-v0) 

I  ;  - - - —  q  d0-.  d0n  ; 

)  J  v  -v  )  cos  4  H1  1 

Sx(a)  1  ° 


1  "W 

1-e  _ 

(v  -v.  )  cos  4 
x  o  X '  To 


q±  d0x  dQ^^  ; 


s1(0) 
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S2( 


l^,(v;  -v0)  1-e-(T/  — V0  ) 


TyT- 

•v0)  cos  % 

(v,  -vj  cos 

(  V-Tj) 

cos  4^ 

l~e 

(▼*«▼,  ) 

(y,-v 

,  )  cos 

(v*.-Y,  )  cos  4e 

(v0-vi)  COS 

*1 

ql  q2  d^l  d6l  d^2  d®2  ’ 

(9a) 


ql  q2  d^l  d9l  d^2  d02  5 


(9b) 


K2(a»0)  ■ 

,a  (fff  J.-'<V'‘VJ 

f  II'/  (v,-v,;  cos  <|'0  **  cos  ^ 

V  2  /  JJjj  - 

S  (a,0) 


(vQ-Vi)  cos  ^ 


ql  ^  d^l  d0l  d^2  d®2  ; 


(10a) 


N~(0,0) 


/V  r2ax,  (v<=-v,  )  cos  4a  (vt-v,  )  cos  <J^ 


s2(0,0) 


(vl-v£)  cos  4> 


qx  q2  d01  d6x  d02  d0£  .  (10b) 


-v 


One  notes  that  e  appears  as  a  factor  in  all  the  transmission  probabilities j 
e  is,  of  course,  the  probability  that  a  photon  is  transnitted  through  the 

slab  without  absorption  or  scattering. 


In  heavy  materials,  the  probability  of  photoelectric  absorption  increases 
very  rapidly  as  the  photon  energy  decreases* 


As  a  consequence,  the 


-  9  - 
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transmission  probability  for  a  path  which  includes  two  or  more  backward 
scatterings  will  be  unimportant  compared  with  the  transmission  probability 
for  a  path  with  the  same  number  of  collisions  all  of  which  are  in  the 
forward  direction;  e.g,  5  (a, a).  In  terns  of  the  integration 

spaces,  the  subspace  corresponding  to  ■  ,n  ■  afl  ■  a  is  the  important 

one  that  need  be  considered  of  all  the  types  Nq*  This  particular  sequence 
of  probabilities  has  the  form 


Nn(a*a»  •  •  *^) 


®n^a,a****a^ 


Vvvi--''.)  ZI  vAAi  • 


on 

where 


f  (v  ) 

ox  o- 


_  W  ~  W  . 

(v  -V,  )  COS 

o  1  o 

fl^v2,Vo^  “  VV.)  • 

*2^V2*Vl,Vo^  "  cos  4^ 

Vw)  -  t2(v2,  V  v0) 
13'V3,V2*vl,vo'  (V2  -v^) cos  f2 


fl(Vo) 


f2(V2'Vl'Vo> 


The  expected  energy,  En  say,  which  is  carried  by  a  photon  trans¬ 
mitted  or  reflected  is  obtained  if  the  integrand  in  Eq.  (1)  is  multiplied 


by  o  ,  where 
J  n* 
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so  that 


E 


n 


n 


-hcsk 


pk+l 


ds,  d0.  ,  d0. 


k+1  k+1 


(12) 


In  Eq.  (12) 

^  . 

pk  “  1  -  (1-cos  ek+]>k 

Every  formula  developed  for  Nn  has  its  counterpart  for  En. 

The  fomulas  summarized  in  Eqs.  (11)  and  (12)  do  not  appear  at  all 
amenable  to  analytic  evaluation.  Numerical  integration,  on  the  other  hand, 
cannot  be  carried  out  for  large  values  of  n  in  a  reasonable  length  of  time. 

It  would  be  of  great  interest,  however,  to  have  a  quantitative  notion  of  how 
many  collisions  contribute  significantly  to  the  number  of  photons  and  the 
energy  transmitted  through  a  large  thickness  of  a  heavy  material.  A  numerical 
calculation  was  therefore  carried  out  for  a  slab  of  uranium  20  cm  thick  upon 
which  photons  are  incident  normally  (<iQ  =  0)  with  energy  aQ  -  1°*  The 
results  through  the  third  scattered  beam  are  given  in  Table  II.  It  is 
estimated  that  these  values  may  be  in  error  by  as  much  as  10*.  The  error 
is  greatest,  of  course,  for  N3  and  E3  and  for  these  the  error  may  be 

pessimistically  estimated  to  be  as  large  as  20*. 

These  results  may  be  used  to  obtain  transmission  values  through  a  slab 

of  lead  for  photons  incident  under  the  same  conditions:  ^  =  0,  aQ  =  10. 

It  may  be  noted  (of.  Table  I  )  that  the  total  absolution  coefficient,  fi(<0, 
for  uranium  is  larger  than  that  for  lead  in  the  energy  range  from  1  to  10  me 
by  an  almost  constant  factor  of  approximately  1.75.  Consequently,  the  integrals 
of  the  types  (7),  (8),  (9),  and  (10)  evaluated  for  20  on  of  uranium  should 
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have  about  the  same  values  for  1«75  x  20cm  «*  35cm  of  lead  except  for  the 
factor  (vr  a/2)  .  The  quantity  in  this  factor  which  varies  in  going  from 
uranium  to  lead  is  /  a,  and  the  ratio  of  1)  a  for  lead  to  its  value  for 
uranium  is  1.09.  In  this  way,  one  gets  the  transmission  values  for  lead 
shown  in  Table  II.* 

These  results  point  to  the  interesting  conclusion  that  the  major 
portion  of  the  photons  and  of  the  energy  trananitted  through  a  slab  of 

heavy  material  of  the  order  of  20  mean  free  paths  thick  is  contained  in 

> 

the  contributions  from  the  first  four  or  live  scattered  beams.  If  an 

u 

extrapolation  is  made  in, a  reasonable  geometric  fashion  to  higher  values 
of  n,  one  gets  for  the  total  energy  and  total  number  of  photons  transmitted 
through  a  slab  of  uranium 


Et  ~  5*5  Eo  ,  • 

(13a) 

N.  r-  10.8  N  , 

O' 

(13b) 

with  the  transmitted  photons  having  the  average  energy 

Et  _  2 

_  5  me  . 

(13c) 

Corresponding  values  for  the  slab  of  lead  are 

Et  ^  6.4Eo', 

(14a) 

Nt  ur  12.7  No  , 

(14b) 

Et  _  ,  2 

n7  -  5  mc  * 

(14c) 

This  extrapolation  is  not  so  daring  as  might  appear.  One  would  expect 

that  once  ER(a)  and  In(a)  begin  to  fall,  they  would  decrease  ever  more 

rapidly  as  n  increases.  This  inference  is  strengthened  by  a  qualitative 

"Note  added  in  proof:  A  direct  calculation  for  35  cm  slab  of  lead  has  been  made 

„  N-,  N  E,  E 

for  aQ  =  10  with  the  following  results:  =  9 .6,  £  =  3 ,h ,  •J  =  1.8u,  -2  =  1 ,8i 

0  ‘^o  ^0  E0 

The  agreement  of  these  numbers  with  the  Pb  values  in  Table  TI  i s  verv  pond. 
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examination  of  the  integrals  for  Ej}  and  In  which  indicates  for  large  a 
a  trend  for  a  behavior  of  the  following  type 


^  (ka) 
*0  •  " 

InM 


n.' 


(ca)n 


n.1 


(15a) 


(15b) 


where  k  and  c  do  not  vary  strongly  with  n.  The  first  four  values  of 

En/Eo  given  in  table  II  agree  fairly  well  with  Eq.  (l$a)  with  the  value 

k  =  ,087cm  the  value  c  =  .12cm”1  in  Eq.  (15b)  gives  a  rough  fit  to  the 

tabulated  values  for  Nn/N0.  Corresponding  values  for  the  slab  of  lead  are  k=.051icm-1j 

c  =  .075cm”1.  It  should  be  noted  that  Eqs.  (l£a)  and  (15b)  are  of  the  form 

which  would  be  obtained  from  a  one-dimensional  calculation  of  the  transmission 

probabilities.  Although  the  values  for  k  and  c  derived  from  a  one-dimensional 

calculation  would  depend  on  the  kind  of  approximation  made,  one  would  not 

expect  that  the  one-dimensional  form  would  underestimate  the  general  trend 


of  the  contributions  to  the  transmission  from  large  values  of  n.  One  may, 
therefore,  argue  that  one  gets  rapid  convergence  in  ^n  and  In,  and  that  an 
extrapolation  made  with  k  and  c  chosen  for  the  best  fit  over  the  first  few 
values  of  n  which  go  beyond  the  maxima  in  ^  and  In,  would  be  reasonably 
safe.  Such  an  extrapolation,  which  uses  Eqs.  (15)  to  get  the  contributions 
for  n  beyond  those  calculated,  gives  essentially  the  same  results  for  E-^ 
and  N-t  as  the  geometric  extrapolation  mentioned  above. 

No  effort  was  made  to  improve  the  accuracy  of  the  computed  values  of 
2n  and  In  presented  in  Table  II  because  of  the  difficulty  involved  in  this 
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straightforward  approach  to  the  calculation.  The  functions  involved  have  a 
behavior  unsuited  to  numerical  methods  and  furthermore  the  number  of  numeri¬ 
cal  quadratures  which  must  be  performed,  for  example,  in  the  case  n  ■  3 
prohibits  the  use  of  more  than  a  few  points  in  each  integration.  The  present 
results  were  obtained  by  using  from  three  to  five  points  for  each  integration. 
So  few  points  would  indicate  a  poorer  accuracy  than  stated  if  one  were  to 
assume  the  use  of  an  unadorned  Simpson's  rule.  The  worst  function  behavior, 
however,  is  found  with  the  functions  p^(0^)  and  qk(©k)  which  are  all  of  one 
class  so  that  they  may  be  considered  together.  Further,  the  functions  of 

0,  and  <U  arising  from  f  (v  v  .  ...  v  )  are  also  more  or  less  of  a  class; 
k  TK  n  n,  n-x,  ,  o 

of  these,  the  greatest  difficulty  for  the  numerical  integrations  comes  from 
a  variation  in  which  the  function  rises  to  a  maximum  and  then  falls  rapidly 
in  roughly  exponential  manner. 

A  typical  integral  which  must  be  evaluated  is  of  the  fom 

I  «V*k(Vdek  •  (16) 

6 

Now  if  the  function  f  (©k)  is  approximated  by  a  polynomial  cQ  ♦  c1©k  +  ... 

+  cn  ©k  which  passes  through  n  points  of  f  and  if  the  integrals 

6k 

J  Q.k  (0^)  d®jc  »  j  ®  0,  1,  2,  ...,  n, 

Gk 

are  evaluated,  then  an  approximate  value  of  (16)  is  obtained  in  which  the 
number  of  points  required  depends  on  the  behavior  of  f  (©k)  and  very  little 
on  that  of  qk  (©k)»  Such  a  procedure  is  useful  since  the  functions  qk 
(and  pk)  are  amenable  to  the  required  analysis.  A  complete  study  of  a  few 
typical  integrals  of  the  type  (16)  showed  that  three  points  can  be  used  to 
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approximate  f  (0^)  with  sufficient  accuracy  to  give  an  error  of  about 

5  percent  in  the  integrated  value*  It  is  essentially  such  a  procedure  which 

the 

was  followed  with  some  care  and  control  to  give/ results  of  Table  II  in  a 
reasonable  time* 

III.  Total  Trftnwnlgsion  as  Successive  Transmission  through  a  Series  of 
Thin  Slabs 

It  is  clear  from  the  proceeding  results  that  the  transmission  through 
a  slab  one  or  two  mean  free  paths  thick  should  not  require  the  evaluation  of 
scattered  contributions  beyond  the  two  collision  beam.  Indeed,  the  twice 
scattered  photons  would  then  constitute  only  a  small  portion  of  the  total 
transmission.  The  calculation  of  the  unscattered  and  singly  scattered  trans¬ 
mission  beams  is  a  simple  matter,  and  the  calculation  of  the  doubly  scattered 
transmission  beam  is  also  not  difficult  if  high  accuracy  in  its  numerical  value 
is  not  required.  With  a  practicable  method  of  evaluating  the  transmission 
through  a  thin  slab,  one  might  propose  to  obtain  the  transmission  through 
a  thick  slab  by  considering  it  to  b®  composed  of  a  number  of  thin  slabs.  This 
procedure  may  readily  be  examined  in  a  formal  manner. . 

Let  0  be  the  operator  which  transforms  a  distribution  of  photons 
incident  upon  one  face  of  a  slab  into  the  distribution  of  photons  leaving 
the  second  face,  and  let  S  be  the  operator  which  transforms  the  incident 
distribution  into  the  distribution  leaving  the  incident  face  of  the  slab. 

It  is  plausible,  in  view  of  the  physical  situation,  to  assume  for  any  dis¬ 
tributions  £  and  r  that  operational  relations 


0  (£  *  r)  m  0£  *  Or, 
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0£  +  5£  ■  (0  +  S)/, 

0  (S/)  -OS/, 

are  valid.  Consider  two  slabs  in  contact  (cf«  Fig.  2)  and  let  be  the 
distribution  entering  the  left  face  of  the  first  slab,  r^  the  distribution 
entering  the  right  face,  and  and  r|  the  distributions  leaving  the  left  and 
right  faces  respectively.  Let  /2,  r2,  /£,  ri,  be  the  corresponding  distribu¬ 
tions  for  the  second  slab.  The  two  slabs,  for  the  general  discussion,  need  not 
be  identical  either  in  material  or  thickness  so  that  the  first  slab  has  operators 

0, ,  S,  and  the  second  3lab  has  corresponding  operators  0  ,  S„.  Then  the 
1*1  A  * 


following  relations  holdi 


1 


°l  h  *  3i  ri 


9 


=1  h  *  °i  ri 


°2  *2  '  S2  r2 


9 


9 


^2  & 2  +  ^2  r2 


» 


and,  since  the  distributions  which  leave  one  inner  face  must  enter  the  other. 


» 


These  six  equations  reduce  to 

ri“°l/l+Slrl  * 

£  ~S1£1*  01  rx  , 

(17) 

r2  “  °2  rl  +  S2  r2  » 
rl  =  S2  rl  +  °2  r2  * 

Substitution  of  the  fourth  of  these  equations  into  the  first  leads  to 

rl  *  ®i  ^1  +  ®l  r2  +  ®1  ^2  rl  * 
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and  this  relation  when  applied  to  itself  gives 

r [  -  (0X  ♦  31  S2  01)  ^  ♦  (31  02  *  S1  S2S1  02)  r2  ♦  S2)2  r^j 

repetition  of  this  process  k  times  yields  the  relation 

“  f  01  ♦  S2  0^  ♦  (S^S2)^DjH%«.  «♦  ( S^S2)  oj 


h  °2  4  S1  S2  Si  °2  +(S1S2)  %*•••♦  <S1S2>  Si  °2  r2 


vk*l  * 

r. 


If  one  substitutes  the  first  of  Eqs.  (17)  into  the  fourth,  and 
iterates  in  a  similar  way,  one  finds 

rl  m[S2  °1  4  S2  S1  S2  °1  +  ^S2S1^2  S2  °1  +”'+  ^S2Sl^k  S2  °1 1  ll 


(IB) 


°2  4  S2  S1  °2  4— 4  <S2S1>  °2 

♦  (s^)*41  r±  . 


(19) 


The  quantities  (SjS,,)^  r,  and  (32^)  rl  rePr®36nt  a  distribution  of 
photons  which  have  been  reflected  back  and  forth  between/ slabs  k+1  times. 

For  k  sufficiently  large,  these  terms  may  be  neglected.  If  Eqs.  (18) 

and  (19)  so  simplified  are  substituted  into  the  second  and  third  of  Eqs.  (17), 

one  finds 

£1  "  !S1  4  °1  S2  °1  4  °1  S2  S1  S2  °1  4  °1  (S2S1)2  S2  0],  +  ** 

U  T  2  1  (20) 

1°1  °2  +  °1  S2  S1  °2  +  °1  (S2S1)  °2  4  —J  r2  * 

r2  ”{°2  °1  4  °2  S1  S2  °1  4  °2  (S1S2)2  °1  4  —  ]  £1 

r  2  1  (21> 

4  iS2  4  °2  S1  °2  +  °2  ^1  S2  S1  °2  4  °2  ^S1S2^  S1  °2  4  **• J  r2  * 
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Eqs.  (20)  and  (21)  give  the  two  distributions  emerging  from  the  two  outer 

of 

faces  of  the  combined  slabs  in  terms/ the  two  entering  distributions.  Eaeh 

sum  of  operator  products  in  the  brackets  of  these  equations  represents  an 

operator  (S  or  0)  for  the  composite  slab  in  terms  ,of  the  operators  for  the 

two  constituent  slabs.  It  is  clear  that  if  a  third  slab  is  brought  into 

the  systan,  the  foregoing  process  could  be  repeated  to  obtain  the  operators 

for  the  crenbination  of  three  slabs*  and  this  process  may  be  extended  for  any 

number  of  component  slabs.  In  any  such  extension,  the  operator  coefficient 

of  will  either  be  the  transformation  of  this  distribution  incident  on  the 

left  face  of  the  composite  slab  into  the  distribution  transmitted  out  of 

the  last  face  on  the  right,  or  it  will  be  the  transformation  which  reflects 

this  distribution  back  from  the  incident  face.  The  behavior  in  a  general 

case  of  any  number  of  component  slabs  is  illustrated  well  enough  by  Eqs,  (20) 

the 

and  (21).  The  coefficient  of  in  Eq.  (21)  consists  of/following  operators: 
0^0^  which  represents  the  operator  for  direct  transmission  through  the  two 
slabs;  C>2  S2  0^  which  represents  transmission  through  the  first  slab,  a 
reflection  from  the  second  slab,  a  reflection  from  the  first  slab,  and 
finally  a  transmission  through  the  second  slab;  etc.  As  would  be  expected, 
this  operator  sum  represents  all  possible  successions  of  transmissions  and 
reflections  which  end  finally  in  transmission.  The  other  operator  brackets 
have  similar  interpretations. 

The  number  of  photons  in  the  distributions  0 2  0^  0 2  5^  S2  0^  .. 

must  be  decreasing.  Indeed,  for  a  material  with  an  appreciable  absorption 
cross  section,  one  would  expect  that  the  important  contribution  to  the 
transmission  would  be  accurately  represented  by  0 0^.  With  this  assumption, 
the  transmission  through  k  identical  slabs  is  given  by  the  simple  operator  0^ 
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This  operational  process  is  a  promising  approach  and  has  the  advantage 
of  giving  details  about  the  photon  distribution  at  a  succession  of  thick¬ 
nesses  of  material.  In  this  way  a  rather  complete  history  of  the  attenuated 
beam  is  obtained. 

The  formulation  in  Section  II  gives  the  necessary  basis  for  the 
construction  of  the  operator  0.  The  distribution  of  photons  incident 
upon  a  slab  may  be  specified  in  terns  of  the  initial  energy,  ag,  and  the 
angle  of  incidence,  The  variables  3  ■  l/ag  and  y  -  cos  are  somewhat 
more  convenient  so  that  the  incident  frequency  distribution  will  be  I  (3,y). 
Then,  for  the  transmitted  distribution,  one  has 

i'  (p'.yVoKm)- 

Now  0  may  be  considered  to  be  the  sum  of  the  operators  0Q,  0^,  0^,  ... 
where  0k  transforms  I(p,y)  into  1^(3', r'),  and  1^  (3',y')  is  the  frequency 
function  for  those  phdtons  which  are  transmitted  with  exactly  k  collisions. 

The  first  of  these  operators,  0Q,  is  found  at  once.  Those  photons  in  the 
incident  distribution  elanent  I  (3,y)  dpdy,  which  are  transmitted  without 
collision,  have  unaltered  enerjy  and  direction  of  travel  and  are  diminished 
in  number  by  the  factor  exp  [-ap(3)/yj,  where  (i(p)  is  the  total  cross  section 
of  the  material  for  photons  of  energy  3  and  a  is  the  slab  thickness.  Hence, 

0Q  is  defined  exactly  by  the  relation 

^  (p\r')  -  exp  i-a+i(fi)/Yj  I(3»y). 

The  operator  0^  may  be  constructed  as  follows.  The  number  of  photons 
of  initial  energy  l/p  and.  incident  on  a  slab  of  thickness  a  with  the  angle 
cos  ^y  which  are  transmitted  after  exactly  one  collision  with  energy  and 
angle  within  the  intervals  ( p * ,3 '  +  dp*  )  and  (y*,y'  ♦  dy')  is 
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dNx  (a) 


V 


r  a  e  -e  _•  9  (8,  0) 

'('Y'/P') 


9  '  U 

y(v-v)  *  5 


1  1 

dy  dp  , 


The  quantities  v  and  v  in  this  expression  are  the  vq  and  respectively 

of  Section  II.  It  follows  that  the  number  of  transmitted  photons  in  the 

11  1  1  1  1 

intervals  (p  ,P  +dp  )  and  (j  ,y  +dy  ),  which  initially  were  in  the  incident 
element  I  (p,y)  dpdy,  is  I(p,y)  dfJdy  •  dN^  (a).  Thus, 


5 


'f 

J  I(i3»Y)y^v.-v) 


-V  -V' 

e  -e 


q’  drd" 


Vr2a 

2 


i(P#y) 


£!=£!' J4  g'l^ci^g'^*2)  p2+(gt3-2p y2)  p;b»2 

v'-v  p73  y(  1-cos ^  0^— y * 2  cos  ^y'y-y^)  V* 


dydp. 


(22) 

where  cos  ■  1  ♦  p-p  .  It  is  evident  that  the  operator  0^  is  a  multi¬ 
plication  by  a  function  followed  by  an  integration  over  the  proper  space 
in  the  variables  p,  y. 

The  operators  O^O^,  etc.  can  be  defined  in  the  same  way.  There 
remains,  however,  the  problem  of  performing  the  indicated  transformations 
with  sufficient  accuracy  and  speed.  The  essential  features  of  one  approach 
to  this  problem  will  be  illustrated  with  the  operator  0^.  The  incident 
distribution  l(p,y)  is  represented  in  the  range  of  interest  or  significance 
by  a  set  of  discrete  points,  i  •  k  in  number,  taken  from  i  values  of  p: 

Pi  <  P2  <  •••<  and  k  values  of  y:  yi  <  y2  <  ...  <  y^.  These  values  of 
I(Psyt)  are  arranged  in  a  column  matrix  // 1 s^/.;  .  Then  I'fp^y*)  as 
determined  by  Eq.  (22)  can  be  evaluated  at  a  given  point  (p  ,y  )  by  performing 
i  numerical  integrations  over  y(j  successively  fixed  at 
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followed  by  a  numerical  integration  over  0.  This  procedure  can  consist  of 
multiplying  each  discrete  value  of  I(0,y)  in  the  column  matrix  by  a  suitable 
coefficient  and  accumulating  the  products.  It  is,  thus,  readily  suggested  that 
one  fora  a  square  matrix  II  0^11  of  order  i  •  k  in  which  the  elements  in  the 
first,  second,  third,  etc.,  rows  are  respectively  the  coefficients  for  ob¬ 
taining  1^  (P^Y^*  l|  (P1# Y2)»  etc*  One  h®8#  consequently,  a 

matrix  representation  for  the  transformation  0^.  A  matrix  representation 
||0  ||  for  Og  can  also  be  found,  and  its  derivation  is  entirely  similar 
to  the  above.  The  operator  0^,  which  is  merely  a  multiplication  by  a  func¬ 
tion  with  no  integrations,  is  represented  by  a  diagonal  matrix.  For  a  slab 
which  is  sufficiently  thin  so  that  the  higher  operators  are  negligible, 
one  has  the  total  matrix 

li  0  II  -  lloji  *  ii  oji  *  ilojl.. 

such  that  the  total  transmitted  distribution  I*  is  given  by 

1 1  I  h  ■  II  0  II  X  //  I  l!  • 

If  an  accurate  calculation  requires  matrices  of  high  order,  the 
feasibility  of  this  procedure  is  brought  into  question.  The  functions 
involved  are  not  well  behaved  which  tends  to  mate  the  order  of  the  matrices 
large.  It  may  be  noted,  for  example,  that  the  integrand  in  Eq.  (22)  becomes 
infinite  at  the  boundary  of  the  region  of  integration.  This  region  is 
defined  by  the  simultaneous  inequalities: 

e'-  2  <  p  <  p'  , 

Y  cos  ©^  -  i  1— y1^  <  Y  <  Y’  cos  ©,  +  1  1-Y'2  )  1-cos^©^  , 

Y  >  0. 

The  area  so  defined  is  clearly  not  rectangular  so  that  there  is  the  difficul¬ 
ty  of  ending  the  numerical  integration  with  respect  to  y  (p  fixed)  on  the 
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boundary  when  no  discrete  point  falls  on  the  boundary.  These  difficulties, 
however,  are  not  insuperable  and  may  be  so  treated  as  to  make  the  method 
feasible. 

The  transmission  of  four  different  distributions  incident  upon  a  slab 
of  uranium  was  calculated  by  this  iterative  method.  The  results  are  presented 
in  Tables  III,  IV,  V,  VI,  and  VII#  Case  A  in  these  tables  represents  the 
first  attempt  with  this  procedure  and  is  largely  an  exploratory  calculation. 
The  incident,  or  initial,  distribution  is  constant  in  the  reciprocal  energy 
range  0.1  <  p  <  0.15  and  is  also  constant  over  the  range  of  oosine  of 
incident  angle  0.8  <  y  <  l.Oj  the  incident  distribution  is  zero  elsewhere# 

A  net  of  49  points  was  used  in  the  computation  and  their  values  are  given 
by  g  ■  0.1,  0.125,  0.15,  0.20,  0.25,  0.35,  0,45  and  y  »  0,  0,2,  0.4,  0.6, 

0.8,  0.9,  1.0.  As  the  calculations  proceeded,  it  became  clear  that  the 
accuracy  could  be  improved  by  a  different  arrangement  of  points  without 
great  increase  in  their  number.  The  improvement  consists  in  taking  a  closer 
spacing  of  y  -  values  for  small  p  and  a  larger  spacing  of  y  -  values  for 
large  p.  Since  the  distributions  which  are  operated  on  vary  rapidly  for 
small  p  and  are  slowly  varying  for  large  p,  it  is  clear  that  this  non- 
uniform  net  is  advantageous.  The  improved  net  was  incorporated  in  the 
matrix  operator  for  the  incident  distributions  which  are  labelled  Cases  I,  ' 

II,  and  III.  For  each  of  these,  the  incident  beam  is  nono-energetic, 

P  -  0.1,  and  the  distributions  over  y  are  respectively  1,  5(y-0.8), 

25  (y-0.8)  in  the  interval  0.8  <  y  <  1  and  zero  elsewhere.  This  choice 
of  angular  distributions  was  motivated  by  several  considerations.  Some 
distribution  of  the  incident  beam  over  a  finite  range  of  y  is  indicated 
for  the  practical  reason  that  a  mono-angular  entering  beam  adds  greatly  to 
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the  computational  difficulties.  On  the  other  hand,  an  incident  distribution 
with  appreciable  magnitude  at  or  near  y  *•  0  would  tend  to  increase  the 
importance  of  the  matrix  11 3  II  which  is  neglected.  The  incident  distribution 
of  Case  I  is  precisely  the  distribution  from  an  isotropic  distant  point 
source  out  of  which  radiation  in  a  cone  of  half-angle  36°  52'  is  accepted. 
The  successive  distributions  of  Case  I,  II,  and  III  may  be  regarded  as  the 
first  three  elements  of  an  infinite  sequence  of  incident  distributions 
which  converges  to  the  mono-energetic  and  mono-angular  beam  (y-1,  o»10) 
considered  in  Section  II. 

The  thickness  of  the  slab  element;  i.e.,  the  "thin"  slab  associated 

with//o  II ,  has  always  been  held  at  2cm.  This  thickness  is  less  than  2 

mean  free  paths,  and  reduces  the  contribution  of  the  twice  scattered  beam 
a 

to /few  percent  of  the  total.  The  ratio  of  the  total  number  of  photons 
transmitted  to  the  number  Incident  is  given  in  Table  III  for  thicknesses  of 
uranium  from  2  to  20cm  in  2cm  intervals.  The  expected  energy  transmitted 
is  similarly  tabulated  in  Table  IV.  It  is  clear  that,  for  the  sequence  of 
Cases  I,  II,  and  III,  the  mean  angle  of  the  initial  distribution  with  the 
normal  to  the  slab  decreases,  and  therefore  the  transmission  should  increase. 
One  may  examine  the  question  as  to  how  rapidly  this  sequence  of  cases  is 
approaching  the  mono-angular  transmission  values.  For  p  ■  0.1,  y  =  1.0  and 

a  total  slab  thickness  of  20cm  of  uranium,  one  has  for  the  ratio  of  the 

number  transmitted  without  scatter  to  the  number  incident.  N  /N,  «■ 

'  o'  i 

exp  ^  -^Qa  !  =  2.89  x  10  ,  This  ratio  and  Eqs.  (13)  give  for  the  mono- 

angular  distribution 

^t  —7  —7 

~  -  3.12  x  10  -  1.59  x  10  '  . 


-  23  - 


p-155 

9  June  1950 


These  values  compared  with  those  of  the  sequences  in  Tables  III 
and  IV  show  slow  convergence  to  the  mono-angular  values*  This  same  slow 
convergence  is  found  in  the  unscattered  transmission  which  can  be  calculated 
simply  and  accurately*  For  20cm  of  uranium,  the  sequence  of  mono— energetic 
cases  considered  here  give  for  the  unscattered  transmitted  number,  Nu,  and 
the  unscattered  transmitted  energy,  Eu,  the  successive  values 
N  E 

"  Ej  "  7,36  x  10’9*  1,14  x  1Cr&l  1*40  *  ID"®- 

The  limit  value  2.89  x.  lCT8  indicates  the  convergence  of  the  unscattered 
transmission  is  fully  as  slow  as  the  convergence  of  the  total  transmission. 

Tables  V  and  VI  summarise  the  values  for  the  ratio  of  the  total  trans¬ 
mission  to  the  unscattered  transmission.  These  ratios  eacpreas  the  trans¬ 
mission  in  a  convenient  manner.  The  u  nscattered  transmission  N  and  E 

u  u 

(or  Nq  and  for  the  mono— energetic,  mono— angular  initial  distributions ) 

are  obtained  by  an  elementary  calculation  which  is  subject  to  no  computational 

inaccuracy;  further,  Nu  and  Eu  give  a  rough  order  of  magnitude  of  the 

transmission.  One  may  compare  the  behavior  of  the  sequence  of  values  of 

w  VEu for  the  cases  I,  II,  and  III  with  the  mono— angular  values  of 

Section  II.  The  values  in  Tables  V  and  VI  Btand  In  the  proper  relationship 

with  each  other  and  with  the  limit.  Again,  the  slow  convergence  of  the 

sequence  of  angular  distributions  toward  the  values  for  the  mono-angular 

case  is  evident.  It  is  worthy  of  note  that  the  effects  of  an  initial 

it  is 

angular  distribution,  even  though/fairly  sharp,  are  very  evident  in  the 
transmission  values,  and  these  effects  persist  as  the  angular  distribution 
is  sharpened.  Another  such  effect  is  to  be  noted  in  the  behavior  of  the 
average  energy  of  the  photons  transmitted.  These  values  are  tabulated  in 


—  2 4  — 
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Table  VII*  For  20cm,  these  average  photon  energies  for  the  sequence  are  In 

o 

proper  relationship  with  eaoh  other  and  with  the  limit  value  of  5  me  for 
the  mono-angular  case* 

The  distributions  for  Chse  3n  of  the  scattered  photons  merging  from  slabs 
of  various  thicknesses  are  shown  in  Fig,  3*  For  purposes  of  comparison, 
these  distributions  have  been  divided  by  the  total  number  of  photons  trans¬ 
mitted  for  the  given  thickness.  Therefore,  the  area  under  each  distribution 
curve  gives  the  percentage  of  the  tdal  transmission  which  comes  from  photons 
which  suffered  at  least  one  collision.  For  slab  thicknesses  of  2,  6,  12,  and 
20<an,  these  percentages  are  20.8,  49*0,  73*1*  and  88.4  respectively. 

The  results  for  uranium  may  be  adjusted  to  give  estimates  for  trans¬ 
mission  through  lead  (Tables  VIII  and  IX).  The  method  follows  closely  the 

one  used  in  Section  II,  It  is  clear  from  the  discussion  given  in  that 
section  that  one  may  write 


CD  00 


where  1^  is  the  appropriate  integral.  The  ratio  of  the  absorption  coefficient 
for  uranium  to  that  for  lead,  it  will  be  remembered,  is  roughly  1.75  in 
the  interval  l.Q  <  a  <  10.0.  Therefore,  1^  for_a  cm  of  uranium  has  approxi— 
mately  the  same  value  for  1,75^®  of  lead.  However,  because  of  the  differ¬ 
ences  in  the  value  of  V  for  lead  and  uranium,  one  must  correct  the  k'th 
tern  by  multiplying  it  by  the  factor  (1.09)  .  In  Section  II,  where  the 
first  few  terms  in  this  sum  were  evaluated  explicitly,  this  correction  was 
straight  forward,  but  here  one  is  given  only  the  sum.  Suppose,  however, 


25 
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one  assumes  that  for  a.  an  of  uranium  N^/?^  -  K/kl,  as  Eqs.  (15)  suggests. 

Then,  if  one  selects  K  so  that 

co 

H  Kk/kl  -  eK  -  N  /N  t 
k-0  *  U 

one  has  for  1.75a  cm  of  lead  W.^AJ ,  -  e^*^k.  The  assumption  basic  to  this 

easy  transformation  of  results  is  probably  not,  in  view  of  its  purpose, 

since 

too  seriously  in  error /  the  transformed  result  is  not  very  sensitive  to 
the  manner  in  which  N./N  is  distributed  over  the  terms  of  the  series.  The 

u  U 

greatest  divergence  from  the  assumption  occurs  for  thin  slabs,  but  for  these 
the  first  tern  which  needs  none  of  this  correction,  gives  most  of  the  total 
transmission. 

As  a  general  remark  concerning  the  tabulated  values,  the  use  of 
three  significant  figures  mif£it  be  taken  as  indicative  of  over-optimism 
regarding  the  accuracy  of  the  computations  because  of  their  extent  and 
complexity.  However,  it  is  estimated  that  a  numerical  result  from  one 
matrix  iteration  is  accurate  to  2%  or  better.  Since  the  error  may  be 
cumulative,  it  might  be  implied  that  10  iterations,  which  carry  the  trans¬ 
mission  through  20cm  of  uranium,  could  have  an  error  of  20%*  Such  an 
accumulation  of  errors  cannot  be  dismissed  only  by  saying  that  it  is  im¬ 
probable.  However,  the  results  of  the  various  calculations  support  each 
other  with  notable  consistency,  and  this  behavior  has  lead  the  writers  to 
the  view  that  the  results  and  the  procedure  have  very  acceptable  accuracy. 
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TABLE  I*  These  values  for  the  absorption  coefficients 
for  lead  and  uranium  were  used  in  the  calculations,  a  -  photon 
energy  in  units  of  me  {  ■  total  absorption  coefficient 

for  uranium;  »  total  absorption  coefficient  for  lead. 
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TABLE  II.  Transmission  values  for  photons  incident 

p 

normally  with  energy  10  me  (aQ  *  10,  “  0).  The  values 

for  uranium  apply  to  a  plane  slab  of  thickness  a  -  20  cm] 
the  values  for  lead  apply  to  a  slab  of  thickness  a  -  35  cm. 
N  is  the  number  of  photons  transmitted  with  exactly  n 
scatterings.  En  iB  the  expected  energy  transmitted  with 
exactly  n  scatterings. 


Transmission  values  for  U 

a  “10,  “0,  a=2Q  on 

0  0 

Transmission  values  for  Pb 

a  “10,  i/  “0,  a»35  cm 

0  0 

n 

E  /E 
n'  0 

0 

1.0 

1.0 

1 

2.6 

1.7 

2.8 

1.85 

2 

3.0 

1.5 

3.6 

!•<* 

3 

2.2 

.88 

2.85 

1.15 
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TABLE  III.  The  values  tabulated  are  the  ratios  of  the 
total  number  of  photons  transmitted,  Nt'  to  the  number  incident, 
IL  ,  for  various  thicknesses,  _a,  of  plane  slabs  of  uranium. 

The  four  cases.  A,  I,  II,  III,  for  which  is  tabulated, 

represent  various  incident  photon  distributions  in  energy  and 
angle  and  are  specified  in  the  text. 


a  (cm) 

A 

I 

II 

Ill 

0 

1 

1 

1 

1 

2 

2.36  x  10"1 

2.22  x  10"1 

2.34  x  10"1 

2.42  x  10""1 

4 

5.20  x  10”2 

4.68  x  10-2 

5.21  x  10-2 

5.50  x  10~2 

6 

1.12  x  10“2 

9.80  x  10”3 

1.14  x  10-2 

1.23  x  10“2 

8 

2.40  x  10**3 

2.05  x  10”3 

2.46  x  10“3 

2.69  x  10-3 

10 

5.07  x  10"4 

4.29  x  10”4 

5.28  x  10~4 

5.85  x  10”4 

12 

1.07  x  lO*"4 

8.97  x  10-5 

1.13  x  10*4 

1.26  x  10"4 

14 

2.24  x  10~5 

1.87  x  10”5 

2.40  x  10~5 

2.71  x  10'5 

16 

4.70  x  10”6 

3.91  x  10”6 

5.08  x  10”6 

5.79  x  10"6 

16 

9.83  x  10-7 

6.14  x  10”7 

1.07  x  10”6 

-6 

1.23  x  10“° 

! 

!  20 

i  1 

_7 

2.06  x  10  ! 

1.70  x  10"7 

—7  ! 

2.26  X  10  ' 

i 

— 

2.61  x  10“7  j 
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TABLE  IV.  The  values  tabulated  are  the  ratios  of  the 
total  expected  energy  trananitted,  to  the  incident  energy, 
E^,  for  the  conditions  of  Table  III. 


a  (an) 


0 


2.05  x  10“x 

1.85  x  10  A 

1.97  x  10  * 

2.03  x  10 

4.14  x  1CT2 

3.45  x  10~2 

3.86  x  10~2 

4.09  x  10’ 

8.34  x  lo”3 

6.50  x  10”3 

7.63  x  10"3 

8.25  x  10' 

1.69  x  10“3 

1.24  x  10"3 

1.52  x  10”3 

1.67  x  10* 

3.42  x  10"4 

2.42  x  10**4 

3.04  x  10*’4 

3.40  x  10 

6.95  x  10“5 

4.74  x  10“5 

6.10  x  10**5 

6.92  x  10 

1.42  x  lo*5 

9.38  x  10”6 

1.23  x  10“5 

1.41  x  10 

-6 

2.89  x  10 

1.87  x  10"6 

2.50  x  Hf6 

2.68  x  10 

_7 

5.93  xlO  ' 

3.75  x  10-7 

5.08  x  10"7 

5.91  x  10 

1.22  x  10~7 

7.56  x  10“8 

1.04  x  10“7 

1.21  x  10 

p-155 

9  June  1950 

TABLE  V.  The  values  tabulated  ape  the  ratios  of  the 
total  number  of  photons  transmitted,  N,  ,  to  the  number  trans— 
mitted  without  scattering,  Nu,  for  the  conditions  of  Table  III. 
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TABLE  VI.  The  values  tabulated  are  the  ratios  of  the 
total  energy  transmitted,  E^,  to  the  energy  transmitted  with¬ 
out  scattering,  E  ,  for  the  conditions  of  Table  III. 


a  (cm) 

A 

1 

II 

hi 

0 

1.0 

1.0 

1.0 

1.0 

2 

1.31 

1.28 

1.27 

1.26 

4 

1.65 

1.61 

1.59 

1.58 

6 

2.05 

2.03 

1.98 

1.96 

8 

2.53 

2.57 

2.47 

2.43 

10 

3.08 

3.25 

3.08 

3.01 

12 

3.72 

4.11 

3.84 

3.73 

14 

4.45 

5.18 

4.78 

4.61 

16 

5.24 

6.53 

5.94 

5.69 

18 

6.16 

8.20 

7.36 

7.02 

20 

7.18  : 

10.3 

9.12 

: 

8.63 
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TABLE  VII.  The  mean  energy  in  the  transmitted  beam, 
W  is  tabulated  for  the  conditions  of  Table  III, 


a  (cm) 

A 

_| 

II 

III 

0 

8.10 

10.0 

10.0 

10.0 

2 

7.05 

8.36 

8.39 

8.39 

4 

6.45 

7.36 

7.42 

7.43 

6 

6.03 

6.63 

6.72 

6.73 

8 

5.71 

6.07 

6.18 

6.21 

10 

5.46 

5.63 

5.75 

5.80 

12 

5.27 

5.28 

5.42 

5.47 

14 

5.12 

5.01 

5.15 

5.20 

16 

4.99 

4.79 

4.92 

4.98 

18 

4.88 

4.61 

4.74 

4.80  j 

20 

4.80 

4.46 

1  4.59 

j 

4.64 
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TABLE  VIII.  The  values  tabulated  are  the  ratios  of 
the  total  number  of  photons  transmitted,  N^,  to  the  number 
transmitted  without  scattering,  Ny,  for  various  thicknesses, 
a,  of  plane  slabs  of  lead.  The  four  cases,  A,  I,  II,  III 
represent  the  same  incident,  distributions  as  were  used  to 
obtain  the  results  of  Table  III. 


a  (cm) 

A 

I 

II 

III 

0 

1.0 

1.0 

1.0 

1.0 

3.5 

1.54 

1.59 

1.56 

1.55 

7.0 

2.19 

2.34 

2.29 

2.26 

10.5 

2.98 

3.39 

3.24 

3.19 

14.0 

3.96 

4.8 

4.5 

4.4 

17.5 

5.1 

6.7 

6.2 

6.0 

21.0 

6.5 

9.3 

8.4 

1 

8.1 

24.5 

8.1 

12.7 

11.3 

10.7 

28.0 

10.0 

17.1 

14.9 

14.1 

31.5 

12.2 

22.6 

19.6 

18.4  ! 

35.0 

14.6 

30.2 

25.8 

24.0 

j 
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TABLE  IX.  The  values  tabulated  are  the  ratios  of  the 
total  expected  energy,  E^,  to  the  energy  transmitted  without 
scattering,  Eu,  for  the  conditions  of  Table  VIII. 


a  (cm) 

A 

I 

II 

III 

1.0 

1.0 

1.0 

1.0 

3.5 

1.31 

1.30 

1.29 

1.72 

1.68 

1.66 

1.64 

2.18 

2.16 

2.10 

2.08 

2.74 

2.79 

2.67 

2.63 

17.5 

3.40 

3.60 

3.40 

3.31 

21.0 

4.2 

4.6 

4.3 

4.2 

24.5 

5.1 

6.0 

5.5 

5.3 

28.0 

6.0 

7.7 

6.9 

6.6 

31.5 

7.2 

9.8 

8.8 

8.3 

35.0 

8.5 

12.6 

11.0 

10.4 

A  schematic  represented 
transmitted  photon  distri 
and  2  and  for  the  comb 


Figure  3  ' 

The  ratio  of  scattered  transmitted  photons  Ns(a)  to  the  total  transmitted  photons  Nt  is 
plotted  as  ordinate  against  the  photon  energy  a  in  units  of  me2  for  the. incident  distri¬ 
bution  of  Case  III.  The  curves  are  given  for  various  thicknesses,  a  ,  of  plane  slabs 
of  uranium,  and  these  curves  are  normalized  so  that  the  area  under  each  curve  equals 
the  fraction  of  the  total  transmitted  number  which  consists  of  scattered  photons. 

These  fractions  are  indicated  as  percentages  on  each  curve. 


